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Magnetic ground state of pyrochlore oxides close to metal-insulator boundary 

probed by muon spin rotation 
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Magnetism of ruthernium pyrochlore oxides 742RU2O7 (A = Hg, Cd, Ca), whose electronic prop- 
erties within a localized ion picture are characterized by non-degenerate t2 g orbitals (Ru 5+ , Ad 3 ) and 
thereby subject to geometrical frustration, has been investigated by muon spin rotation/relaxation 
(/tiSR) technique. The A cation (mostly divalent) was varied to examine the effect of covalency 
(Hg > Cd > Ca) on their electronic property. In a sample with A = Hg that exhibits a clear 
metal-insulator (MI) transition below ~ 100 K (which is associated with a weak structural tran- 
sition), a nearly commensurate magnetic order is observed to develop in accordance with the MI 
transition. Meanwhile, in the case of A — Cd where the MI transition is suppressed to the level 
of small anomaly in the resistivity, the local field distribution probed by muon indicates emergence 
of a certain magnetic inhomogeneity below ~ 30 K. Moreover, in Ca2Ru207 that remains metallic, 
we find a highly inhomogeneous local magnetism below ~25 K that comes from randomly oriented 
Ru moments and thus described as a "frozen spin liquid" state. The systematic trend of increasing 
randomness and itinerant character with decreasing covalency suggests close relationship between 
these two characters. As a reference for the effect of orbital degeneracy and associated Jahn- Teller 
instability, we examine a tetravalent ruthernium pyrochlore, TI2RU2O7 (Ru 4+ , 4d 4 ), where the result 
of /iSR indicates a non-magnetic ground state that is consistent with the formation of the Haldane 
chains suggested by neutron diffraction experiment. 

PACS numbers: 75.25,+z, 75.50.Lk, 76.75. +i 



I. INTRODUCTION 

A class of transition metal oxides with a chemical 
composition A2-B2O7 crystalizes into so-called cubic py- 
rochlore structure (Fd3m), and it can be regarded as 
two interpenetrating networks of corner-shared tetrahe- 
dra consisting of A 2 0' and B2OQ sublatticesJ- Because 
of high crystal symmetry, presence of antiferromagnetic 
correlation or charge disproportion between the metal 
ions on the pyrochlore lattice induces a geometrical frus- 
tration and associated macroscopic degrees of degeneracy 
for the spin/charge states in their ionic limit. This causes 
a variety of interesting phenomena such as a "spin glass" 
phase in i? 2 Mo 2 7 with R = Y, Tb and Dy (RefsHi) 
and the "spin-ice" phase in i? 2 Ti 2 C>7 with R = Dy and 
Ho (Refsill). 

It is naturally expected that geometrical frustration 
(or correlation) would also lead to nontrivial effects on 
itinerant electrons, because they can be a local entity 
in a certain time scale (~ h/st, where £t is the rele- 
vant energy band width). The unusual behavior of or- 
dinary and anomalous Hall coefficients in the Mo py- 
rochlores with i?=Nd, Sm and Gd might be a manifes- 
tation of such effect£i£ Above all, however, the "heavy 
Fermion" behavior in a metallic vanadium spinel LiV 2 04 
(in which the V ions form a pyrochlore sublattice) is one 
of the most interesting phenomena among conducting 
pyrochlores/spinels4 

In the previous studies on LiV 2 04, we demonstrated by 



muon spin rotation/relaxation (^SR) and muon Knight 
shift measurements that, unlike what is anticipated in 
the conventional heavy Fermion metals, local vanadium 
moments in LiV 2 04 are not quenched by the Kondo 
mechanism, and that they give rise to a highly inho- 
mogeneous and dynamically fluctuating magnetic ground 
state that seems to be in close association with the ge- 
ometrical frustration^*^ It would be interesting to note 
that a strikingly similar situation was reported to occur 
in an intermetallic Laves phase (C15-type) compound, 
Yi_ a; Sc 2 ;Mn 2 , where it is inferred from earlier /iSR study 
that inhomogeneous local magnetism persists at low tem- 
peratures in a sample that exhibits a heavy Fermion be- 
havior induced by suppression of structural phase transi- 
tion upon substitution of Y by ScJ^ Thus, these results 
demonstrate a pressing need to develop a new concept of 
quasipartcle state that should represent the anomalous 
conducting state of electrons under geometrical (three- 
body) correlation. 

From the experimental view point, it would be helpful 
for the understanding of anomalous metallic properties 
to establish the generic character of such electronic states 
by extending /iSR study to other electrically conducting 
pyrochlores/spinels. To this end, spinels and pyrochlores 
are good candidates for materials survey, although they 
exhibit a common tendency of structural phase transition 
at lower temperatures that lifts the orbital degeneracy 
and relieve the frustration. 

Here, we focus on a family of ruthenium and osmium 
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FIG. 1: Single ion energy manifold of B site ions in 
^2 + ^2 + 07 with B=Ru, Os. Three d-electrons are distributed 
over different t^g orbitals (d xy , d yz , d zx ), so that no orbital 
degeneracy is expected. 



pyrochlores, A2-B2O7 (B=Ru, Os) with nominally diva- 
lent cations (i.e., Al + Bp~0 7 ). As shown in Fig. [TJ the 
B site ions have 4d 3 (5d 3 ) configuration with the d or- 
bits split into e g doublet and t2 S triplet states (due to 
the crystal fields of BOq octahedrons), where three elec- 
trons are distributed to different t2g orbits (xy, yz, zx) 
to satisfy both Pauli exclusion principle and Hund cou- 
pling. This means that there is no orbital degeneracy 
that would lead to the instability against the Jahn- Teller 
distortion and associated structural phase transition to 
reduce the crystal symmetry, and therefore one may be 
able to observe the effect of geometrical frustration at 
lower temperatures. Moreover, considering that the sup- 
pression of structural transition is favorable for prevent- 
ing the system from MI transition, they would provide 
a stage for studying the effect of geometrical correlation 
on itinerant electrons. Interestingly, the actual ground 
state of these pyrochlores suggested from bulk electronic 
properties turns out to be considerably dependent on the 
choice of A site ions. As summarized in Table U some 
exhibit metal-insulator (MI) transition that does not nec- 
essarily accompany structural phase transition. We have 
shown by recent /iSR study that an osmium pyrochlore, 
Cd2 0s20j, which remains in a cubic structure and falls 
into semi-insulating state below 225 K, exhibits a pecu- 
liar magnetic ground state probably related with geomet- 
rical frustration.^- A similar but more inhomogeneous 
electronic state is suggested by a preliminary /xSR study 
for the case of Hg20s207 which remains metallic down 
to 4 K (RefOJ). 

In this paper, we report on the magnetic ground state 
of ruthenium pyrochlore oxides (A2RU2O7) that are sit- 
uated near the MI transition. For the nominally pen- 
tavalent ruthenides (id 3 ), the A cation (~ A 2+ ) is var- 
ied among Hg, Cd, and Ca to examine the effect of co- 
valency (Hg > Cd > Ca) on their electronic property. 
It is shown in samples with A = Hg and Cd that a 
nearly commensurate magnetic order is observed to de- 
velop in accordance with the MI transition (below ^100 
K, which is associated with weak structural transition). 
However, the local field distribution at the muon site in- 



dicates emergence of a certain magnetic inhomogeneity in 
Cd2Ru 2 07 below ^30 K, suggesting that the spin state 
retains significant amount of entropy at lower temper- 
atures. Meanwhile in Ca2Ru207 that remains metal- 
lic, we find a highly inhomogeneous local magnetism 
below ^25 K that is described by a Gaussian distri- 
bution for the local fields at the muon site character- 
istic those from randomly oriented magnetic moments 
situated on a regular lattice. The observed trends of 
greater randomness in the magnetism with increasing 
itinerant character of (i-electrons is common to that ob- 
served in osmium pyrochlores, and thereby it suggests 
a close link between the degree of magnetic randomness 
and that of delocalization. As a reference for the effect 
of orbital degeneracy and associated Jahn- Teller insta- 
bility, we also examine a nearly tetravalent ruthernium 
pyrochlore, TI2RU2O7 (Ru 4+ , 4<i 4 ), where the result of 
jiiSR indicates a non-magnetic ground state that is con- 
sistent with the formation of the Haldane chains sug- 
gested by nuclear magnetic resonance (NMR) 14 and neu- 
tron diffraction experiments ! 15 ! 16 



II. EXPERIMENT 

Polycrystalline samples of A2RU2O7 (A =Hg, Cd, Ca, 
and Tl) used in this study were prepared by a solid-state 
reaction of AO and Ru02 together with oxidizer at 950 
C° under a high pressure of 4 GPa. The details of sample 
preparation were described elsewhere^ The bulk prop- 
erties of these samples were investigated by variety of 
methods including magnetic susceptibility, electric resis- 
tivity, and specific heat measurements. Structural anal- 
ysis was made using powder X-ray diffraction technique 
(combined in part with synchrotron radiation). Result 
of these analyses are summarized in Table HI In order to 
obtain sufficient quantity of samples for good signal to 
noise ratio in /iSR measuremnts, samples with A =Cd, 





SPT (T s ) 


MIT (Tmi) 


X-drop (T m ,T g ) 


Hg2Ru2C>7 

Cd 2 Ru2 07 
Ca2Ru207 


~ 107 K 
No? 
No 


~107K 
No? K 

No 


~ 107 K 

~ 100 K 

SG for T g < 25 K 


T12RU207 


- 120 K 


~ 120 K 


- 120 K 


Hg20s2C>7 

Cd 2 Os 2 07 


No? 
No? 


No 
~230 K 


~ 80 K" 
~ 230 K±2 



TABLE I: Electronic ground state of Ru and Os pyrochlores, 
where SPT refers to structural phase transition suggested 
by X-ray diffraction data, MIT to metal-insulator transi- 
tion suggested by resistivity data. Magnetic susceptibil- 
ity (x) in Hg2Ru207 and Cd2Ru207 exhibits a drop below 
Tmi, suggesting development of antiferromagnetic spin corre- 
lation, whereas a spin glass-like (SG) hysteresis is observed in 
Ca2Ru207. (The temperature dependence of x are shown in 
the following sections together with the ^iSR result.) 



Ca, and Tl were prepared in two batches under the same 
condition and mixed for exposure to muons after char- 
acterization by bulk property measurements. Their sus- 
ceptibility data exhibited a slight variance of transition 
temperatures (i.e., 2-3 K for A =Cd and Ca, ~ 5 K for 
A=T1). 

Among these compounds, Hg2Ru20y is characterized 
by a sharp reduction of magnetic susceptibility (%) be- 
low T m ~107 K with a small hysteresis. This occurs in 
conjunction with increase of resistivity (p) by nearly one 
order of magnitude and small lattice expansion accom- 
panying splitting of X-ray diffraction peaks attributed 
to simple cubic unit cell.il These obserations led to an 
earlier speculation that spin-singlet ground state might 
be realized in association with the structural phase tran- 
sition. It is interesting to note at this stage that the 
tetravalent TI2RU2O7 shares many features in common 
with Hg2Rri2 07: it also exhibits sharp reduction of \ a t 
T m ~ 120 K that is associated with a jump of p and 
structural phase transition . 15 ! 16 

On the other hand, Cd2Ru20y does not show any clear 
sign of structural change in the lattice parameters around 
T m , although it exhibits reduction of x below T m ^100 
K that is similar to Hg2Ru20j. Resistivity exhibits a 
rather peculiar temperature dependence, where it is char- 
acterized by an anomalous kink around T m and divergent 
increase at lower temperatures that suggests insulating 
ground state. Similar features have been reported in an 
earlier literature.— 

Compared with these two compounds, Ca2Ru2 0y is 
more distinct in the sense that it remains metallic at 
least above ^2 K; although p{T) shows slight enhance- 
ment at lower temperatures, it does not obey the rela- 
tion p(T) oc ex.p(E a /kT) expected for thermally acti- 
vated carriers over a band gap E a . More interestingly, \ 
exhibits a cusp around T m ~ 25 K that is typically found 
in spin glass systems. The present sample, which was 
synthesized by solid-state reaction under high pressure, 
exhibits these features in common with that previously 
obtained by hydrothermal method. 19 

Conventional /iSR measurement was performed using 
the Lampf spectrometer installed on the M15/M20 beam- 
lines at TRIUMF, Canada. During measurements under 
zero external field (ZF), residual magnetic field at the 
sample position was reduced below 10 -6 T while the ini- 
tial muon spin direction was parallel to the muon beam 
direction [P M (0) || z]. For longitudinal field (LF) mea- 
surement, a magnetic field was applied parallel to P (U (0). 
Time-dependent muon polarization [G z (t) = z ■ P«(t)] 
was monitored by measuring decay-positron asymmetry 
along the z-axis, 



A(t) = A Q G z (t) 



N + (t)-aN_(t) 



N+(t) + aN-(ty 

where, Aq is the average asymmetry, 

N ± (t)=N ± (0)e~ t ^{l±A ± G z (t)} 



(1) 



is the positron event rate for the detector placed in the 
forward (+) or backward (— ) position relative to the sam- 
ple, r M is the muon decay lifetime (= 2.198 x 10~ 6 s), 
A± is the decay positron asymmetry for respective de- 
tector (A± ~ Aq), and a is the instrumental asymme- 
try [a = N+(0)/N-(0) ~ 1 in usual condition]. Trans- 
verse field (TF) condition was realized by rotating the 
initial muon polarization so that Pfi{0) || x, where the 
asymmetry was monitored along the x-axis using signals 
from an appropriate pair of detectors for N± (£) to obtain 
G x {t) = x ■ P„{t). 

Because of the limited quantity of samples (0.1-0.2 g) 
obtained by synthesis under high pressure, a special de- 
vice ( "muon veto counter" ) was installed to reduce back- 
grounds due to positron events associated with muons 
that missed the sample. All the measurements under a 
magnetic field were made by cooling the sample to the 
target temperature after the field equilibrated. 



III. RESULT 
A. HgaRuaOy 

Fig. [5] shows some examples of ZF-/zSR time spectra 
over a temperature range from 2 K to 210 K, where one 
can observe abrupt emergence of oscillatory signal around 
T m ~ 109 K without enhanced depolarization due to crit- 
ical slowing down of spin fluctuation in approaching the 
transition temperature. As is also evident in the Fast 
Fourier transform (FFT) of these spectra shown in Fig. [3l 
they indicate step-like development of internal field just 
below T m associated with a long-range magnetic order. 
These features are consistent with suggestion from the 
hysteresis of x(T) that the magnetic transition [which is 
not obvious solely from x(T)] is of the first-order, taking 
place in cooporation with the structural phase transition. 

The FFT spectra are characterized by two bands of 
dominant frequency lines around 10 MHz and 30 MHz at 
2 K, suggesting possibilities of internal field distribution 
coming from i) a commensurate spin density wave probed 
at a unique muon site or ii) an antiferromagnetic state 
with several crystallographically non-equivalent muon 
sites. Considering these possibilities, we analyzed the 
time spectra using the following model function valid for 
the ordered phase of polycrystalline samples, 



A(t) =j Y^A, 



1 



-Al* 



G^ 1 (i,A) e - A ' t cos(27r/ 4 i 



3 3 

+A b , (2) 

which, in the paramagnetic phase, is reduced to 

A(t) = (A -A b )Gf T (t,A) + A b , (3) 

where G* T (t, A) is the Kubo-Toyabe relaxation 
function- to describe the slow Gaussian depolarization 
due to random local fields exerted from nuclear magnetic 
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FIG. 2: (Color online) Time-dependent /i-e decay asymmetry 
[A(t) oc muon polarization] observed in Hg2Ru207 at various 
temperatures under zero external field, where the full polar- 
ization corresponds to ~ 0.21. Development of spontaneous 
internal magnetic fields is inferred from the sinusoidal oscil- 
lation of A(t). 
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FIG. 3: (Color online) Fast Fourier transform of the ^SR 
spectra shown in Fig. [2] Each spectrum is shifted vertically 



by 2.5 x 1CT 4 
clarity. 



(or 5 x 1(T 4 between 105 K and 107 K) for 



moments (with A ~ 10 _1 MHz being the linewidth in 
the quasistatic limit), Ai is the partial asymmetry, Al is 
the longitudinal relaxation rate (which turned out to be 
~ for Hg2Ru2C>7), Ai is the transverse relaxation rate, 
fi is the muon spin precession frequency (= ^7^2?$ with 
7 M = 2tt x 135.53 MHz/T and B, being the local field at 
the muon site) , 4> is the initial phase of precession, and 
Af, is the background coming from muons that missed 
the sample E^i + Ab = A ). When the dynamical 
fluctuation of local fields is negligible, we have 



GT(t,A) 



§♦§<■ 



i 



A 2 i 2 )exp(--A 2 £ 



2.2\ 



(4) 



exp(-AV), (At < 1) 



for which the original density distribution is given by 



P KT (B a ,A) = 



7a* 



exp 



7 2 B 2 
' 2A 2 



(a = x,y,z). (5) 



For the longitudinal component in Eq. (J2J), one must 
keep in mind that the one third of the implanted muons 
are subjected to a local field parallel to the initial muon 
polarization in the magnetically ordered phase of poly- 
crystalline samples, so that the depolarization due to nu- 
clear random local fields are suppressed (A <C 7/j-Bi)- We 
adopted four components (n = 4) which turned out to be 
minimal for the satisfactory description of the entire spec- 
tra, although the spectra turned out to be dominated by 
two of these (i = 1 and 3, see below). 

As shown in Fig. [4j the step-like development of the 
internal magnetic field below T m is inferred from the tem- 
perature dependence of fi and of the fractional yields for 
the magnetically ordered phase (Ai/Ao). In particular, 



the total fraction shown in Fig. H[b) indicates that nearly 
100% of the entire crystal falls into fully ordered state 
just below T m : a small background (A b /A ~ 5%) is at- 
tributed to potassium chloride that remains as a residue 
product of oxidizing agent (KCIO3). The depolarization 
rate of the precession signals is relatively large for the 
major two components (A; ~ 27r/i). Moreover, while it 
is mostly independent of temperature above ~ 50 K, it 
exhibits a tendency of increase with decreasing temper- 
ature at lower temperatures (see Ai) that is also visible 
in the FFT spectra. As is seen below, this tendency be- 
comes much significant in the case of Cd2Ru20y. The 
small value of the longitudinal relaxation (Al) indicates 
that the local field is static in the time window of obser- 
vation (< 10 _1 MHz), and that the large value of A; is 
primarily due to the random distribution of local fields 
at the muon site (not to the dynamical fluctuation). 

The Gaussian linewidth, A, is determined by local 
configuration of the nuclear magnetic moments ( Ru, 
101 Ru, 199 Hg, and 201 Hg in this case) which are nearest 
neighboring to muons, and therefore it provides useful in- 
formation for identfying muon sites. The observed value 
[A ~ 0.04(1) MHz] has been compared with calculated 
dipole sum, 



rSBiVi 



1/2 



(6) 



for several possible muon sites, where A 3 - is the dipole 
tensor, 



1 



3ai/3j 



<W (a,/3 = x,y,z), (7) 
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FIG. 4: (Color online) Temperature dependence of bulk sus- 
ceptibility (a) and fit parameters obtained by analysis using 
Eq. p). namely (b) fractional yield for the ordered phase 
defined by Ai/Ao, (c) central frequency (ft), and (d) depolar- 
ization rate (Ai) in Hg2Ru207. 




0.2 2 4 6 
Time ((is) 



FIG. 5: (Color online) Time-dependent /i-e decay asymmetry 
observed in Cd2Ru2C>7 at various temperatures under zero 
external field. 



is 90°, they suggest a magnetic structure of multiple-q 
(7^ 0) in which Ru spins are antiparallel between those 
mutually at the diagonal sites of Ru-hexamers on the 
Kagomc lattice (perpendicular to the [111] axis). Our 
simulation of internal field at the presumed muon site 
with further assumption of antiferromagnetic spin struc- 
ture for triangular lattice layers reproduces the two-band 
structure of muon precession frequency with the local 
fields ~ 0.2 and ~ 0.4 T/hb- The comparison of these 
values with the experimental result yields estimated Ru 
moment size of ~ 0.5(1) hb, indicating the considerable 
shrinkage of the moment size form 3/is. 



and the sum is made for the contribution of the j'-th 
nuclear magnetic moments Ij located at Tj = (xj,yj,Zj) 
from a given muon site. As a result, it turned out that 
an interstitial position near the center of Ru-cornered 
tetrahedra is the most probable site (calculated value of 
A = 0.04 MHz). 

Based on this muon site assignment, one can further 
proceed to examine its consistency with the proposed 
magnetic structure by comparing the magnetic fields ex- 
erted from the Ru 5+ electronic moments 



B, 



loo 



= \J2^ 



Mil 



(8) 



with the internal field B t experimentally observed (~ 0.2 
T for i = 1 and —0.1 T for i =2, 3), where fij is the mag- 
netic moment of Ru ions. Recent nuclear magnetic reso- 
nance (NMR) study using Hg 199 and enriched 99 Ru nuclei 
suggests a rather complicated magnetic structure consist- 
ing of four non-equivalent Ru sites.— Considering that i) 
the internal field is canceled at a quarter of Hg nuclear 
sites, and that ii) the angle between the primary axis 
of the electric field gradient and Ru magnetic moment 



B. Cd 2 Ru 2 7 

The ZF-/iSR spectra in Cd2Ru207 are qualitatively 
similar to those observed in Hg2Ru20y. A signal showing 
spontaneous precession of muon polarization emerges be- 
low T m ~ 100 K, indicating that a magnetically ordered 
state develops below T m . The absence of enhancement in 
the relaxation rate with decreasing temperature toward 
T m is a feature common to Hg2Ru207. It is inferred from 
the Fourier transform of these spectra that they consist 
of two bands of frequency lines, which are most obvious 
around 80 K. However, such a feature, once fully fledged 
over the intermediate temperature, becomes gradually 
blurred with decreasing temperature below ~ 50 K. This 
is clearly observed by comparing the FFT spectrum at 
2 K, where it shows only a single broad band while the 
corresponding spectrum in Hg2Ru207 retains two-band 
structure (although the latter also exhibits a common 
tendency of broadening below ~ 50 K). It is also notice- 
able that the muon precession frequency is considerably 
reduced from -15-28 MHz in Hg 2 Ru 2 7 to -12-19 MHz 
in Cd2Ru207, suggesting that the effective Ru moment 
size is accoordingly reduced. 
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FIG. 6: (Color online) Fast Fourier transform of the ^SR 
spectra shown in Fig. [5] Each spectrum is shifted vertically 
by 2.5 x 10~ 4 (or 7.5 x 1(T 4 between 80 K and 95 K) for 
clarity. 



Curve fits of these spectra in time domain using Eq.Q 
reveals that we need three components, of which one re- 
mains paramagnetic (/3 = 0) down to 2 K. The tem- 
perature dependence of fit parameters is summarized in 
Fig. [71 While the internal field develops sharply below 
T m (as seen in Fig. Wp), the fractional yield of the mag- 
netically ordered state [(^4i + A 2 )/A , Fig. [TJd] exhibits 
more gradual increase with decreasing temperature than 
that observed in Hg 2 Ru207 (Fig. |4pl). Moreover, it ex- 
hibits a decrease below T m i ~ 30-40 K in place of the 
paramagnetic component (A3). The change is also associ- 
ated with enhanced depolarization of the first component 
(i = 1) whose frequency merges into that of the second 
component (i — 2). It is also noticeable that the bulk 
susceptibility exhibits a kink near T m 2 in correspondence 
with the occurrence of strong depolarization. 

The small relaxation rate for the longitudinal depolar- 
ization (Al, not shown in Fig. [7J at 2 K indicates that 
the origin of the enhanced depolarization for the preces- 
sion signal (A2) is mainly due to the enhanced quasistatic 
randomness of the local fields at the muon site. This has 
been also supported by the recovery of the asymptotic 
component in /iSR spectra observed under a longitudi- 
nal magnetic field (LF). However, it might be worthy 
of mention that a broad peak of Al has been observed 
around 40 K, which suggests slow residual fluctuation of 
local Ru moments at the intermediate temperature. 

It is inferred from the value of A [=0.04283(2) MHz 
which is close to the case of Hg2Ru2C>7] that muon sites 
in Cd2Ru207 are common to those suggested for the case 
of Hg 2 Ru2C-7. Then, the observed change of the /uSR 
spectra below T m i is ascribed to that in the magnetic 
structure of Ru moments, which might be associated with 
the suppression of MI transition in Cd2Ru207. Assum- 
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FIG. 7: (Color online) Temperature dependence of (a) bulk 
susceptibility (corresponding to one of the two-batch samples) 
and fit parameters obtained by analysis using Eq. ([2}, namely 
(b) fractional yield for the ordered phase defined by Ai/Ao, (c) 
central frequency (/i, where fz — 0), and (d) depolarization 
rate (Ai) in Cd2Ru207- 



ing that the magnetic structure is common to that of 
Hg2Ru2C>7 for T m <2 < T < T m , the reduction of internal 
field at the muon site is attributed to the shrinkage of 
effective Ru moment size to 0.36(7)/i_b, which is consis- 
tent with enhanced itinerant character of d-electrons in 
Cd 2 Ru 2 C>7. 



C. Ca 2 Ru 2 7 

As mentioned earlier, the bulk property of Ca2Ru207 
is markedly different from the two previous compounds, 
where the magnetic property inferred from susceptibil- 
ity is close to that expected for the classical spin glass 
systems. The ZF-^tSR spectra in this compound also 
exhibits a number of features that are known as those 
characteristic for the spin glass. As found in Fig. |5Ji, 
they exhibit increase of depolarization rate and associ- 
ated change of line shape from Gaussian to exponential- 
like damping with temperature approaching to the spin 
glass transition temperature (T g ~ 25 K), indicating 



the critical slowing down of spin fluctuation towards T g . 
However, the detailed analysis reveals that the depolar- 
ization is not described by a simple exponential function 
expected for the case of unique fluctuation rate. More- 
over, the line shape below T g (which is determined by 
the distribution of quasistaic local magnetic fields) turns 
out to exhibit a distinct feature which is not seen in the 
canonical spin glass like dilute alloys. 

It is well known that the density distribution of in- 
ternal magnetic field (and associated hyperfine coupling) 
in the dilute alloys (CuMn, AuFe, etc., with Mn and Fe 
concentration less than 3%) is described by a Lorentzian 
distribution, 



with that for electronic moments, 



P h {B a ,a)=P±)- 2 

V -k I a 1 



llBl 



(a = x,y,z) (9) 



with a being the factor determining the width of 
distribution ! 22 ' 23 This comes from the strong random- 
ness of RKKY interaction which is effectively described 
by 



s(r) 



cos(kpr) 
(kpr) 3 



(10) 



where fcp is the Fermi momentum, and r is the distance 
between the local magnetic moments. The distance r 
varies randomly in the dilute limit while the interaction 
is long-ranged (decays only by the third power of r), lead- 
ing to freezing of randomly located spins along random 
orientation. The corresponding time evolution of muon 
polarization under zero external field is known to have a 
form 

G L z (t,a) = l + l(l-at)e X p(-rt), (11) 

4 
~ exp(--cri). (crt< 1) 

o 

Apart from the influence of spin fluctuation, the observed 
/iSR spectra may be then simply written as 

A(t) = A Gf T (t,A)G^(t 7 <j) + A b , (12) 

which would represent the line shape at the lowest tem- 
perature (T <S^ T g ). Here, we assume that the hyperfine 
coupling between the electronic spins of Ru atoms and 
nuclear spins are small enough so that the spin dynamics 
of the latter is independent of the former, and that they 
remain quasistatic within the time window of /zSR. This 
allows us to make an approximation that the muon depo- 
larization is described by a product of the Kubo-Toyabe 
function (nuclear-spin part) and exponential decay (elec- 
tronic part) [which is also implicitly assumed for Eq. (|2|)], 
However, as shown in Fig. [5b, the comparison of Eq. (TT2")) 
with the actual spectra observed at 1.9 K indicates that 
this form does not reproduce the data. More interest- 
ingly, the line shape at 1.9 K is excellently described by 
the Kubo-Toyabe function with the line width replaced 



As 



± 7 2 £(|A^| 



1/2 



(13) 



where Ag ^> A as \fij\ 3> \Ij\. The comparison in Fig. |SJd 
is made for a set of spectra with different external longi- 
tudinal fields (Ho) to enhance the asymptotic component 
in Eq. (|lll) [where these curves can be obtained numer- 
ically by replacing B z in Eq. ^ with B z — [j,qH , and 
the same is true for Eqs. (@| and ([3]) with A replaced 
by As]. The response of asymptotic component, particu- 
larly the dip around t ~ 0.02-0.04 /xs and quick recovery 
to the asymptotic value, indicates that the line shape 
is described by the quasistatic local fields that obey a 
Gaussian distribution. 

We also find in the detailed analysis that the observed 
spectra are reasonably reproduced only when we assume 
multiple components with different fluctuation rates of 
the local fields, so that Eq. (TT2|) should be eventually 
replaced with 

n 

A(i)^ G? T (i, A) ^G^i, A s , ^) + A, (14) 

where G^ T (t, Ag,i/) is obtained numerically by solving 
an integral equation for the "strong collision" mode l 20 ' 24 , 

Gr(t,A s ,v) = G™{t,A s )e- vt 

+u f e-Gf(r, A s )Gf(t - r, A s , v)dr, (15) 
Jo 

which is deduced from the assumption that spin-spin cor- 
relation is described by 



m ■ m 
\m\ 2 



(16) 



The best fits are obtained with two components [n = 2 in 
Eq. (|14p] for the spectra below T g with a Gaussian width 
A s = 38.5(7) MHz (as found in Fig. |SJ the asymptotic 
component exhibits a slow depolarization corresponding 
to 2/2 = 2.3(8) x 10 8 s" 1 while v x < 10 5 s" 1 at 1.9 K), 
whereas three components (n = 3) are required to yield 
reasonable fits for these above T g . In such a situation, 
simultaneous fits of multiple spectra obtained under dif- 
ferent longitudinal magnetic fields using a common set 
of parameters provide reliable results. One such set of 
data obtained at 30 K yields result of curve fit shown in 
Fig.El with v x = 1.48(18) x 10 9 s~\ v 2 = 1.72(34) x 10 8 
s _1 , and fs = 0.69(11) x 10 7 s _1 for a common Gaussian 
width As = 19.3(5) MHz, all having comparable partial 
asymmetry. As shown in Fig. 1101 these values of spin 
fluctuation rate spread over three orders of magnitude, 
strongly suggesting that the spin dynamics of Ca 2 Ru 2 07 
cannot be reduced to a simple model like Eq. ([TBI as- 
suming a unique correlation time (r c = 1/v) at a given 
temperature. 
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FIG. 8: (Color online) (a) Time-dependent fi-e decay asym- 
metry observed in Ca2Ru207 at various temperatures under 
zero external field, and (b) the initial part of the time spectra 
at 1.9 K under a longitudinal field of 10 mT, 50 mT and 100 
mT. Solid curves are the best fits using Eq. (|14[) assuming the 
Gaussian distribution for the local fields (leading to the Kubo- 
Toyabe relaxation function), whereas dashed curves are those 
obtained by the Lorentzian distribution (with a — As/\/2). 



0.20 



Ca 2 Ru 2 7 LF-p.SR(30K) 



0.15 



< 



▼ LF=210mT 

a LF=100mT 

• LF= 50mT 

o LF= lOmT 



0.10 A 



0.05 



0.00 




8 10 



Time (u,s) 



FIG. 9: (Color online) (a) Time-dependent fi-e decay asym- 
metry observed in Ca2Ru207 at 30 K under a longitudinal 
field of 10 mT, 50 mT, 100 mT and 210 mT. Solid curves are 
the best fits using Eq. (|14[) assuming three components. 
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It turns out that the nuclear dipolar width in 
Ca 2 Ru 2 07 [A = 0.128(2) MHz, determined from the 
spectra for T 3> T g without much ambiguity] is con- 
siderably greater than those observed in Hg2Ru2 07 and 
Cd2Ru 2 07, suggesting that muon seems to occupy a dif- 
ferent site closer to Ru ions: 43 Ca is the only isotope of Ca 
nuclei with finite nuclear magnetic moments whose nat- 
ural abundance is negligibly small (0.135%) and there- 
fore of no relevance to the change of A. However, con- 
sidering that both nuclear and atomic dipolar moments 
are situated in the same Ru site, the magnitude of As 
can be used to evaluate the effective atomic moment size 
of Ru ions from that of nuclear magnetic moments us- 
ing Eqs. © and JI3J), yielding 0.69(l)/i B /Ru (T < T g ) 
or 0.35(l)/x B /Ru (T > T g ). The latter is in excellent 
agreement with the value estimated from the tempera- 



FIG. 10: (Color online) Temperature dependence of (a) bulk 
susceptibility (corresponding to one of the two-batch samples) 
and (b) Fluctuation rate of the local fields probed by /iSR in 
Ca2Ru207. Three components with different fluctuation rates 
(^1,2,3) are presumed in Eq. ()14|) for T > T g ~ 25 K, whereas 
two of them (^1,2) for T < T g (see text for the detail). 



ture dependence of bulk susceptibility obtained for single- 
crystalline specimen of Ca2Ru2C>7 J£ 

It might be appropriate to point out at this stage that 
the above mentioned result clearly indicates the differ- 
ent origin of randomness in Ca2Ru207 from that in the 
dilute alloy systems. This is readily understood by con- 
sidering the difference in the nature of spin-spin inter- 
action between these two cases. The RKKY interaction 



in the dilute alloys has a specific feature that its sign 
changes randomly with relative distance between the lo- 
cal moments located at random. Meanwhile, in the case 
of pyrochlore compounds, all the magnetic moments are 
located at the regular lattice position, where the random- 
ness is allowed only for the orientation of the moments: 
there is no obvious route to introduce randomness to the 
interaction between B-site ions. The fact that the local 
field is described by a Gaussian distribution in Ca 2 Ru2 07 
also indicates the random orientation of Ru 5+ local mo- 
ments (S — 3/2) as the primary source of the random- 
ness, which might be phenomenologically described as a 
"frozen spin liquid" state. 




4 6 

Time ((is) 



D. TI2R.U2O7 

Unlike other three compounds, the Ru ions in 
TI2RU2O7 are presumed to be tetravalent (4d 4 ) in the 
ionic picture, and therefore one of the Ad electrons has 
orbital degeneracy in the £29 manifold. Such a degen- 
erate state is susceptible to the Jahn- Teller distortion 
that reduces the free energy by breaking the local lattice 
symmetry, generally leading to structural phase transi- 
tion. As a matter of fact, TI2RU2O7 is known to exhibit 
a clear structural transition at T s ~ 120 K in conjunc- 
tion with the metal-insulater (or -semiconductor) tran- 
sition, where the reduction of symmetry from cubic to 
orthorhombic structure is inferred from powder neutron 
diffraction studiesi 15 ' 16 The transition is also associated 
with a sharp drop of magnetic susceptibility (x) a t ^s, 
where the observed small hysteresis suggests that the 
reduction of \ might be driven by the structural tran- 
sition and associated opening of the band gap. These 
bulk properties are quite similar to those observed in 
Hg 2 Ru2 07, leading to a naive speculation that the elec- 
tronic ground state might be common between these two 
compounds. 

However, it has been suggested from Tl-NMR experi- 
ment that the ground state of TI2RU2O7 is non-magnetic, 
and that it has an excitation gap (as inferred from the be- 
havior of the spin-lattice relaxation rate, T^ 1 ).— The re- 
cent inelastic neutron scattering experiment also reports 
occurrence of a well-defined excitation gap (~ 11 meV) 
below T s »i£ The gap is attributed to the formation of one 
dimensional Ru chains (S = 1) along [110] direction that 
are understood as a model system of Haldane chains^ 
Our /zSR result supports the presence of non-magnetic 
ground state, and it provides further evidence for the 
formation of the Haldane chains in this compound. 

Fig. [TT] shows fiSK time spectra under zero external 
field observed at various temperatures. They do not ex- 
hibit any anomaly in passing through T s ~ 120 K, where 
the line shape is mostly like a slow Gaussian depolariza- 
tion expected for the random local fields from nuclear 
magnetic moments. (More specifically, the finite gra- 
dient dA(t)/dt for t — > indicates that an exponential 
damping due to paramagnetic fluctuation overlaps to the 



FIG. 11: (Color online) Time-dependent \i-e decay asymme- 
try observed in TI2RU2O7 at various temperatures. Each 
spectrum is vertically shifted by 0.02 for clarity. Inset: the 
asymmetry observed at 2 K under a longitudinal field of 0, 10 
mT, 50 mT, and 100 mT. Solid curves are the best fits using 
Eq. JTT]). 



Gaussian damping; see below.) Thus, the result clearly 
demonstrates that the new ground state emerging be- 
low T s does not accompany long-range magnetic order, 
as has been suggested by Tl-NMR. Meanwhile, one may 
notice that the line shape exhibits a gradual change from 
Gaussian to exponential damping with decreasing tem- 
perature, and that a small fraction of the asymmetry 
shows rapid depolarization below T* ~10 K. This qua- 
sistatic magnetism can be attributed to staggered mo- 
ments at the open edges of Haldane chains,— where the 
edges may come from various crystalline defects including 
the twin boundaries inherent in the orthorhombic struc- 
ture. Some earlier /iSR study on another Haldane system 
(Y2BaNi05, which happens to have a similar gap energy) 
indicates that such staggered moments induce spin glass- 
like random magnetism. 28 

Considering these features, we made curve fit analysis 
of the spectra using a simplified form of Eq. (|T2l with 
two components, 

2 
A(t) ~ GK T (i,A)^A 4 cxp(-A, ; i)+A 6 , (17) 



A, 



2ct 2 ,^ 



(t = l,2) 



which is a good approximation in the paramagnetic state 
with Vi > Cj. The hyperfine parameters (itj and A) are 
determined by the simultaneous fit of multiple spectra 
obtained at 2 K under different longitudinal fields (Ho), 
yielding a x = 0.49(1) MHz, a 2 = 8.3(3) MHz, and A = 
0.042(7) MHz. The value of A indicates that muon site 
is common to those in Hg2Ru 2 07 and Cd2Ru207. These 
values are fixed for the rest of analysis to extract the 
spin fluctuation rates (i/,-) from curve fits for the ZF-/1SR 
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FIG. 12: (Color online) Temperature dependence of (a) bulk 
susceptibility (corresponding to one of the two-batch sam- 
ples) and fit parameters obtained by analysis using Eq. (|17|) . 
namely (b) fluctuation rate (vt) and (c) corresponding depo- 
larization rate (\i). 



spectra, where the result is summarized in Fig. 1121 The 
second component is observed only below T* with A2 
showing a gradual increase with decreasing temperature 
(~ 0.03 at 2 K). 

Interestingly, the fluctuation rate exhibits a gradual 
decrease as temperature is reduced toward T s , which is in 
good correspondence with the behavior of Tl-NMR T^ 1 
interpreted as a manifestation of reduced spin excitation 
associated with the development of spin-singlet correla- 
tion. The fact that it levels off below T s supports the 
view that /iSR is also monitoring the same phenomenon 
through a different time window (i.e., it is not the artifact 
associated with muon diffusion, for example). 

Further slowing down of fluctuation (v\) is observed 
below T* , where the second component attributed to the 
staggered moments appears with intermediate fluctua- 
tion rate (1/2). It turns out that the emergence of this 
component occurs in parallel with the appearance of a 
large Curie term in the magnetic susceptibility. Thus, 
the seemingly "impurity" -like behavior in \ is explained 
by the edge states associated with the Haldane chains. 



Hg2Ru2C>7 

Cd 2 Ru 2 7 



107 K 
100 K 



Ca 2 Ru 2 7 T g ~ 25 K 



0.5(l)/xs 

0.36(7)a*b 
0.35(l) M s (T>T g ), 
0.69(1) fiB (T < T„) 



TABLE II: Effective Ru moment size in A 2 Ru 2 C>7 estimated 
from the magnitudes of local fields at muon site. 



IV. DISCUSSION 



Correlation between magnetic randomness and 
itinerant character of 4d electrons 



We have shown by the /LtSR studies on a series of Ru 
pyrochlore compounds (A2R.U2O7) that their electronic 
ground state exhibits a systematic change with varying 
A cations. The change in the structural and MI transi- 
tion temperature suggests that the relevant 4d electrons 
become more itinerant in the system with less covalent 
A cation (Hg > Cd > Ca). This is also supported by 
the effective size of local Ru moments summarized in Ta- 
ble II, which shows tendency of reduction along with the 
above mentioned order (except below T g in Ca2Ru20y, 
where the moment seems restored on the suppression of 
spin fluctuation) . While it is not clear how the variation 
in the covalency of ^40 sublattice may affect the Ru0 6 
network, the tig band width (ede) seems to be compa- 
rable to (or greater than) the energy gain brought by 
the local lattice distortion. We might stress here that 
the magnetism in the ground state exhibits stronger ran- 
domness in the compounds showing more itinerant char- 
acter. The absence of structural change upon MI tran- 
sition in Cd2Ru207 might even suggest that the tran- 
sition is of a Mott type (U > £d c , with U being the 
electronic correlation energy), as the ti g band is half- 
filled. The appearance of more random magnetism at 
lower temperatures (below ^30 K) suggests that the 
magnetism just below T m ~ 100 K is not a fully fledged 
Neel state, and that a significant amount of entropy is 
preserved. In this regard, it might be worthy of remark 
that the Dzyaloshinski-Moriya (DM) interaction would 
serve as a direct driving force for the secondary magnetic 
phase with an energy scale of 25-30 K (i.e., ~ ksTms, in 
Hg2Ru207, Cd2Ru207, and ~ fcgTj in Ca2Ru20-7 which 
are close with each other and suggesting a common ori- 
gin). The path of exchange interaction between Ru ions 
indeed lacks inversion symmetry, and thereby Ru ions are 
subject to the DM interaction that tends to rotate the 
Ru moments around the direction of Si x Sj . This may 
be cooperatively serving with the itinerant character of 
d electrons (see below) to induce the randomness in the 
electronic ground state of Ru pyrochlorcs. 

According to the two-dimensional t-J model on a tri- 
angular lattice (in the vicinity of half-filled band, where 
"t"~ Ede is the hopping integral, J is the Heisenberg 
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exchange energy) 29 i 30 , the presence of geometrical frus- 
tration leads to the competition in the ground state be- 
tween the Nagaoka's ferromagnetism (which is rigorously 
proven for J = with one hole) and the Neel state 
for t < 0, leading to the large degeneracy and associ- 
ated entropy near T = 0. The strong randomness in 
the magnetism observed in Cd2Ru2 07 and Ca2Ru2 07 
might be understood as consequence of a similar situ- 
ation (with decreasing effective J in place of t) on the 
three-dimensional stage, where the quasistatic character 
of the local spins might imply a small effective J due to 
the two competing magnetic correlations. 



B. Close similarity between Ca2Ru207 
Lii_ a; Zn a: V204 



and 



It has been known in the earlier studies that substitu- 
tion of Li + with Zn 2+ in L1V2O4 induces spin glass phase 
at low temperatures [for x > 0.05 in Lii_. E Zn a ;V204, 
where the valence state of the vanadium is Y 3 ' 5 ~ x ' 2 in 
the ionic limit] j2ir— This observation comprises the ex- 
perimental basis that the geometrical frustration is at 
work among vanadium ions which also occupy the py- 
rochlore sublattice in the spinel structure. The com- 
pound exhibits metallic conductivity except x > 0.9 
where a structural transition (cubic-to-tetragnal) occurs 
at low temperature (~ 50 K at x — 1). Our detailed ^fSR 
study on Lii-^Zn^A^C^ has shown that the fluctuation 
of vanadium electronic moments (> 10 9 s _1 ) which per- 
sists down to the lowest temperature in pure LiV204, is 
strongly suppressed below 4-5 K for x = 0.05-0.1, and 
that a quasistatic random magnetism shows up below 
T g ~ 10 K for x > 0.2£ The sensitiveness of V spin 
fluctuation to the slight variation of band filling suggests 
that the quater-filled state (3d 1 ' 5 for x — 0) in the vana- 
dium t^g band may be a "quantum critical point" (QCP) 
across the competing ground states. 

Here, we point out that the result of /iSR in Ca2Ru20y 
shows a marked resemblance with that obtained in 
Lii_ a; Zn a ;V204. As shown in Fig. [Tj|J the line shape of 
/iSR spectra observed in the quasistatic spin glass phase 
of Lii_ a; Zn a ;V204 with x = 0.2 is well represented by the 
Kubo-Toyabe relaxation function (where the hyperfine 
parameter As — 60 MHz) with a small fraction (~ 20%) 
of correction term [an exponential damping added to sim- 
ulate Eq. (jlip probably due to the randomness induced 
by Zn substitution for Li], indicating that the density 
distribution of local magnetic field probed by muon is 
mostly represented by a Gaussian distribution^ This 
again suggests that the distribution is determined by the 
random orientation of V moments, and that the origin of 
randomness is common to the case of C^Ru^Oy. Consid- 
ering that the spin glass state is connected to the heavy 
quasiparticle state in LiV204, the electronic ground state 
of Ca2Ru207 might be also understood as that close to 
the QCP realized in L1V2O4. 
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FIG. 13: (Color online) Time-dependent /i-e decay asymme- 
try observed in Lii_ I Zn a; V204 with x — 0.2 at 1.6 K under a 
longitudinal field of 0, 10 mT, 0.1 T, and 0.28 T. Solid curves 
are the best fits using Eq. (|14[1 plus a small fraction of additive 
exponential damping term (see text for detail). 



V. CONCLUSIONS 

We have shown by muon spin rotation study on a se- 
ries of ruthenium pyrochlore oxides that the electronic 
ground state exhibits a systematic variation for differ- 
ent A-cations, where the structural and metal-insulator 
transition temperature (T s ) is reduced with increasing 
covalency of the AO sublattice. The reduction of T s and 
enhanced itinerant character of the compounds is in ac- 
cordance with the degree of randomness in the ground 
state magnetism probed by /iSR, where the randomness 
comes from random orientation of the ruthenium local 
magnetic moments. Considering a theoretical investiga- 
tion using the two-dimensional t-J model on a triangu- 
lar lattice, we point out that the origin of such "spin- 
orientarion glass" state may come from the competition 
between the Nagaoka's ferromagnetic correlation and an- 
tiferromagnetic exchange interaction in the vicinity of 
the Mott-insulating phase (for the half-filled T^ g bands) 
which is in co-operation with the local Dzyaloshinski- 
Moriya interaction. 
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